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Introduction. Hamiltonian Formalism of Analytical Mechanics has been 
systematically used after Poincare especially by people who created Quan- 
tum Mechanics in the 20's. In pure mathematics, formalism of differential 
forms appeared as a by-product of Hamiltonian Theory formalized finally by 
E.Cartan. However, geometrical understanding of many important parts of 
Hamiltonian Formalism has not been elaborated on for a long period. For 
example, general definitions of such fundamental geometrical objects as La- 
grangian Submanifolds in Symplectic (Hamiltonian) linear spaces (and in 
more general nonlinear symplectic manifolds as well) were finally formulated 
only in the 60's. In particular, V.Arnold participated in this (see [jl], |2|). 

In the late 60's, the present author observed that some algebraic version 
of Hamiltonian Formalism for rings with involutions plays a fundamental role 
in many constructions of Differential Topology (see 0). As a reaction to this 
work, I.M.Gelfand observed (and pointed out to me in about 1971) that von 
Neumann's construction of selfadjoint extensions for symmetric operators is 
based in fact on some Lagrangian planes in Hilbert spaces with symplectic 
scalar product. 

The present work is in a sense a natural continuation of these old ob- 
servations. An essential part of its results already was announced by the 
present author in a short note [|J.Let me point out here that we started at 
first to discuss graphs with A.Veselov as a continuation of a series of papers 
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dedicated to factorization and Laplace transformations of the Schrodinger 
operators (see || ||). 

Discrete Schrodinger Operators on Graphs 
Wronskians and Topology 

We shall consider Graph T-i.e., one dimensional simplicial complex such 
that: 

1. Only a finite number of edges Ri (equal to nip) can meet each other 
in any vertex P; 

2. Graph T has no ends, i.e. nip > 1 for every vertex P. Two spaces of 
scalar complex-valued functions will be used: 

The space of complex- valued functions ipp depending on vertices P and the 
space of functions depending on edges R. We do not formulate any global 
restrictions on these functions now. 

We shall consider real self-adjoint (i.e. symmetric, in fact) operators L 
acting on these spaces of functions: 

(Lip) P = bp.p'ijjp', b P:P/ = b p/:P (1) 

P:P' 

(Lip) R = d R :R'ipR>,d R .. R , = d R > :R (2) 

R:R' 

Reality means that all coefficients are real. 

Definition 1 Operator L will be called Finite Type iff for any point P or 
edge R there exists only a finite number of points P' or edges R' such that the 
coefficients above are nonzero. The Operator will be called an operator of Fi- 
nite Order iff this number does not depend on P or R. Order of operator 
is a maximal number of vertices ( edges ) in the minimal paths joining such 
pairs of vertices (edges) for which coefficients are nonzero. For the Second 
Order Operators on the vertices nonzero coefficients can be only bp.p = Vp 
and bp : p> iff P\JP' = dR. In the case of edges nonzero coefficients can be 
only d R . R = V R and d R , R i iff Rf]R' ^ 0. We call coefficients Vp and V R 
Potentials 
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For any simplicial complex K and number k we have a boundary operator 
d from A;-chains into k — 1-chains. We also have a scalar product where 5- 
functions (i.e. functions taking value 1 on one of /c-simplices only) give an 
orthonormal basis in subspace of finite chains (i.e. finite functions from k- 
simplices). Therefore we have an adjoint coboundary operator d*. 

Combinatorial Laplace-Beltrami Operators A& on the spaces of k- 
chains are defined by the general formula 

A k = {d + d*) 2 = dd* + d*d 

For the cases k = and k — 1, Laplace-Beltrami operators on Graph T 
are the second order operators (Schrodinger Operators on the vertices and 
edges) such that 

A :b P:P , = 1,P\JP' = dR,V P = -m P (3) 
A 1 : d R:R , = 1, R H i?V 0, V R = -2 (4) 

Example 1 Let Graph T is a line R 1 , i.e., 1-simplices R and vertices P can 
be numerated naturally by the integers n such that boundary of the edge with 
number n is equal to the union of vertices with numbers n and n — 1 . We 
always have mp = 2 here. There is a natural one-to-one correspondance here 
between edges and vertices with the same numbers such that A = Ai . Both 
operators can be considered as the same operator L acting on the functions 
on the Lattice Z . 



(L ip) n = ^ n _! + ip n+1 , ne Z (5) 
-L Q = A + 2 = A 1 + 2 

Two standard bases of solutions for the equation L ip = \ip can be given 
by the obvious formulas 



^± = a™,a± = l/2(A±v / A^4) (6) 

n ) >~>n 

a + — a_ a + — a_ 

C = 1,C 1 = 0,S = 0,S 1 = 1 

^± = C + a±S (7) 
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More general operators L on the Lattice Z appeared in the discretized 
Theory of Solitons (theory of Toda Lattice) beginning in 1974 (see in the 
survey and encyclopedia articles J7], [|). They have form of the second order 
selfadjoint Schrodinger Operators on the graph T = R 1 above 

(Ll[)) n = Cn-lVW-l + V n if) n + c n ip n+1 (8) 

Solving the most fundamental problems of Spectral Theory for these oper- 
ators for rapidly decreasing coefficients c n — l,V n and for periodic (quasiperi- 
odic) coefficients, the so-called Wronskian for any pair of solutions has been 
used. 

Our goal is to invent and to use a natural analog of this quantity for any 
finite order Schrodinger Operator on the arbitrary Graph T satisfying to the 
conditions above. 

Consider any pair of solutions Lip = AiJj, Lip = \<f) on the Graph T. 

Definition 2 The following quantity will be called Wronskian for any pair 
of solutions. 
For vertices: 

R 

W R (4>, ip) = 6 P:P /(0p^p/ - 4> P >(f)p) (9) 
P\JP' = dR 

For Edges: 

R 

W R ((j), V>) = d R:R'((pRtpR' ~ 1/JrM (10) 
Ri 

Rf)R> = P 

We can see that our formula for the Wronskian is very much the same as 
for the line in the case of vertices, but for the edges it is less obvious. For 
example, it contains summation along the edges R' meeting our edge R in 
one point P only. Second boundary point P' where OR = P U P' does not 
appear in the sum. Therefore, in this case the correctness of our definition 
should be proved. 
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Theorem 1 For any pair of solutions 0, ip of the second order difference 
equation L0 = A0, Lip = \ip, Wronskian is a well-defined 1-chain W on the 
Graph Y (i.e. a complex-valued function of the oriented edges R = PP') 
whose boundary is equal to zero 

dW = (11) 

Therefore our Wronskian belongs to the first homology group H = H° pen (T, Z) 
modulo infinity (if Graph is noncompact) . It defines an H -valued skew sym- 
metric scalar product on the spaces of solutions. 

Proof. In the case of vertices, Wronskian is obviously well-defined as 1- 
chain. From the equality (L(p) P ip P — (Lip) P (p P = we extract immediately 
that 

E W R' = 

P'UP=dR' 

Here the edges are taken with such orientation that they end in P. How- 
ever, this equality means precisely that dW = by definition of the boundary 
operator. 

Consider now the case of edges. Our formula above defines correctly this 
quantity as Wr^p depending on the edge R and its vertex P. Starting from 
the same equality (L(f)) R ip R — (Lip) R (p R = for the pair of solutions, we can 
see that the left-hand part is obviously equal to the sum of two expressions. 
One of them is precisely equal to W Ry p as it was defined above for the case 
of edges. The second one is equal to the analogous sum W Rt p' with point P 
replaced by P' . So we have 

W R , P + W R , P , = (12) 

Therefore we have W R:R = —W Rj p>. Wronskian is well-defined as a 1- 
chain. Consider now the boundary of it. 

(dW) P =Y,W R = (13) 
Pei? 

R'^R R 

R'f)R = P,<j ) ' = 0^,0 = ^ 
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However, in the last sum any fixed pair R, R' appears twice with opposite 
signs. Therefore the total sum is equal to zero. Our theorem is proved. f\ 

Higher Order Operators: For the case of higher order Schrodinger 
operators acting on vertices, we define Wronskian in the same way as for the 
second order operators. For any pair of interacting vertices P, P', we fix one 
simple path 

I PP , = R 1 ,R 2 ,..., R k ,dl = PUP' (14) 

Elementary 1-chain associated with interacting pair PP' is defined as 
before. Full Wronskian is a sum of these expressions along all interacting 
pairs of vertices: 

W = E w h Wj = b P:P/ ((j)p^p - iI)p4p) (15) 

dW = 

As before, it is easy to prove that 

W e H^ en (T,Z) 

For the complex solutions ip P , one may consider Wronskian W(ip,ip) as 
a Quantum Current; its Topological property to be 1-cycle is a Kirchhof 
Law. 

Multidimensional Simplicial Complexes: For the natural classes of 
Schrodinger operators acting on the spaces of fc-chains in Simplicial Com- 
plexes K with k > 1, we can define analogous quantity (Wronskian) as a 
function on the set of all pairs Ws k> s k -i (here we have a fc-simplex and its 
k — 1-face). It is exactly the class of operators for which fc-simplices can 
interact if they have common (k — l)-face. This is an exact definition of this 
class of operators and of the Wronskian: 

{I4)s u = E E H-s^s* (16) 
s k ^\£S k s' k ns k =s k -i 

w Sk , Sk .M^)= t E W^V-W (i7) 

s' k r\S k =s k -i 

tp' = ^S k ntp = ^5 fc 

3 It is strange , but the present author was not able to find this elementary fact in the 
literature (1 asked several experts in the theory of graphs and operators on them). It does 
not surprise me for the edges, but I cannot believe that this fact is new for the case of 
vertices. 
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The following theorem is true. 

Theorem 2 1. The full sum of "Wronskians" along all k — 1- faces of every 
k -simplex is equal to zero; 

2. The full sum of "Wronskians" along all k-simplices with common 
k — 1-face is also equal to zero. 

Proof of this theorem is exactly the same as for k — 1. 

Remark 1 For nontrivial applications of the theorems above we need to con- 
sider such classes of operators that for some A the space of solutions is at 
least 2- dimensional, otherwise our Wronskian will be identically equal to zero. 
There are two natural sources for that for the graphs: 

1. Graph F and operator L admit nontrivial symmetry group. 

2. Graph T has several number of "Tails" at infinity and operator L has 
asymptotically constant coefficients. It is exactly a case for which Scattering 
Theory can be developed. We shall do this in the next paragraph. In par- 
ticular, we shall demonstrate that "Unitarity" of Scattering follows from the 
Topological Property of Wronskians established above. 



Scattering Theory for Graphs with Tails 

and 

Symplectic Geometry 

We consider graphs with k tails here. It means precisely that outside 
of the finite domain this Graph is isomorphic to the union of k positive 
half lines ("Tails") z 1: z 2 ■ ■ . , z k . Let us choose some initial vertices Pj,j = 
1,2, ... ,k in these tails and attribute to them the number of the tail and 
zero. Other vertices in the same tail behind Pj will be naturally numerated 
by the number of the tail and positive integers n G Z + . We attribute to the 
edge R = (n, n — 1) number n, as in paragraph 1 for the line, n — 0, 1, 2, 

Any Graph Y with k tails can be naturally represented in unique way as 
a union of the finite subgraph without ends V C V with some number of 
Trees attached to it in the vertices Qi, I = 1, . . . , s. Obviously we have s < k 
because every tree contains at least one tail (maybe more). 
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Definition 3 We call Graph T' the Basis of our Graph T. Points Qi from 
which Trees grow up in the Graph T will be called Nests. Several trees can 
grow up from one nest. Connected Graph T is Topologically trivial iff its basis 
T' is equal to one point. We call graphs T with k tails Diagrams for which 
Scattering Processes naturally can be defined. 

Consider now a class of second order Schrodinger Operators L on Graph 
T with k tails acting on the vertices or on the edges such that outside of the 
finite domain we have in the tails for n > n and all j = 1, . . . , k: 

-L = -L = Ai + 2,i = 0,1 (18) 

Outside of this finite "domain of interaction" our equation Lip = Xijj has 
solutions in the tails ipf n and Cj n , Sj n for j = 1, . . . , k and n > n , described 
above (see formulas (6,7) in paragraph 1). 

Definition 4 Symplectic space H 2k with basis Ci, Si, . . . , Ck, St and real- 
valued skew scalar product such that 

<! Cj, Cp Sj, Sp >= 
<C j ,Sp>=5 jp ,j,p,= l,...,k (19) 

will be called an Asymptotic Symplectic Space. 

Solutions ipf in the tails can be expressed as complex linear combinations in 
the same basis with natural linear extension of the same scalar product. 

Definition 5 Vector ip ass e H 2k such that there is a continuation of it as 
a solution ip on the whole Graph Y will be called an Asymptotic Vector 
for the operator L; \-dependent linear space T\ of all asymptotic vectors 
ip ass E T\ for the operator L will be called Space of Symplectic Scattering 
Data for L. 

Theorem 3 For any self adjoint real second order Schrodinger Operator L, 
the corresponding linear space of Symplectic Scattering Data is Lagrangian 
subspace in the Asymptotic Symplectic Space for every complex value of A. 
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Proof. We shall demonstrate that this fact is, in fact, topological. It 
follows from Theorem 1 which claims that Wronskians are homological cycles 
in our Graph for every pair of solutions <fi, tp. 

At the same time, we know the following elementary topological proper- 
ties of tails in the connected graphs: 

l.Any piece of tail Zj may appear as a nonzero part of 1-cycle only if all 
edges of this tail belong to this cycle with the same coefficient; 

2. Individual tail Zj never has a continuation to Graph T as a cycle con- 
taining other tails with coefficients equal to zero. The difference between any 
pair of tails be continuated to Graph T as a cycle such that other 

tails do not participate in it. 

For the Wronskian of two solutions on T we have 

W((j), ip) = J2 a i z i + (f^ite) = (20) 

p=k 

= P P (zi - z P ) + {finite) 

p=2 

From this equality we conclude that 

E«i = (21) 
i=i 

At the same time we know that this sum is exactly a skew symmetric 
scalar product of asymptotic vectors: 

j=k 

Y^aj =< <f) ass ,^ ass > (22) 
j=i 

Therefore our theorem is proved. 

Our Schrodinger Operators will be considered now as a selfadjoint opera- 
tors in the Hilbert Spaces of square integrable functions L l 2 (T) for two cases: 
vertices i — and edges % — 1 as before. The subset of 5- functions 5p G L\ 
or 5r G L\ contains functions equal to 1 on some vertex (edge), and zero 
otherwise. 

For trivial reasons we always have a continuous spectrum with 
multiplicity equal to k for |A| < 2. Therefore the results of para- 
graph 1 are very useful here. 
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Sometimes we may have exceptional discrete eigenvalues for |A| < 2 drown 
in the continuous spectrum (see below). A "Normal" discrete spectrum ap- 
pears for |A| > 2 

Definition 6 We call the Scattering Zone of Spectrum the area where 

|A| < 2, A G R. The interval of real |A| > 2 we call the Zone of Normal 
Discreet Spectrum, for which every eigenfunction has exponential decay 
in all tails and is nonzero at least in one tail. Exceptional Discrete Spec- 
trum may appear for any real A. It is such that corresponding eigenf unctions 
are identically equal to zero in all tails. 

In the Scattering Zone we have 

ipf n = expjim^A)} = a±,n > n (23) 
0(A) G R, |A| < 2 

In the Zone of Normal Discrete Spectrum we have 

6(\) G iR,i 2 = -1 

So in the last zone both solutions ip^ are real. One of them has exponen- 
tial decay for large n: 

if)j~ — > oo, n — > oo, j = 1, . . . , k (24) 
ipj — > 0, n — > oo 

Let us denote the open halflines on the real A line by 

I + = (2 < A < oo), 7_ = (-2 > A > -oo) (25) 

Using the result of the previous Theorem we have a pair of mappings 

T+ : /+ - A fc , T~ : /_ -> A k (26) 

defined by the Lagrangian Planes T\ for both cases A G I±. Here is a 
Lagrangian Grassmanian whose points are Lagrangian fc-planes in H 2k . 
There is a canonical codimension 1 cycle (see ||) 

ZdA k 

containing all Lagrangian planes with a nonempty intersection with a set of 
directions J2j KjipJ ■ You may say that this cycle contains all singularities of 
the projection, when you project Lagrangian planes into the linear span of 
all vectors ipf forgetting halfbasis ipj. 
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Corollary 1 Normal discrete eigenvalue appears exactly where the Lagrangian 
Plane T\ crosses the cycle Z . Therefore two "Morse Indices" are defined; they 
characterize topological properties of normal spectra for A > 2 and A < 2: 
they are the intersection indices of the curves T 1 * 1 with the cycle Z defined 
above. 

We can see that decay of any eigenfunction belonging to the Hilbert space 
should be exponential here (or it should be equal to zero in all tails- see 
later). 

Remark 2 / did not clarified yet whether these intersection indices are ex- 
actly equal to the numbers of eigenvalues (i.e. all crossing points for the 
generic case have sign +) or not. 

Let us define now Scattering Matrix for Operator L. Consider the Scat- 
tering Zone |A| < 2 where our basic complex solutions 

ij)f n = exp{±in6 j (X)} 

are complex adjoint to each other. Take them as a basis for the complexified 
Symplectic Space H% k where Lagrangian subspace T A is given by the theorem 
above. Generically this subspace can be interpreted as a graph for the linear 
map from halfbasis ip + = to the halfbasis ip~ = (ipj) 

^ - E *jp1>7 (27) 
j 

^+->S(A)ty") 

This corresponds to the choice of basis in the form 

ej=^ + + E s ^7 e7 A (28) 
j 

Corollary 2 In the Scattering Zone the Scattering Matrix S\ is always uni- 
tary S G Uk and symmetric S f = S 

Proof. For k = 2, the relationship between Lagrangian planes and real uni- 
modular matrices was mentioned as an example in the elementary textbook 
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of Arnold (see ||): From the Lagrangian plane in H A , we come to the stan- 
dard Monodromy Matrix M G SLi2(R) which maps basis C±, S\ into the basis 
G-iiSi- This matrix was always in use in Classical Math Literature for the 
second order Sturm- Liouville operator on the line. Coming to the com- 
plex bases ipf and ^fj=l,2, we get a monodromy matrix in the group SU\ t \ 
which is isomorphic to SL 2 (R). Monodromy matrix is unimodular because of 
wronskian property. Algebraic transform from this to the Scattering Matrix 
has the standard name of Caley Transformation in this case. It was used 
by quantum physicists in the Scattering Theory for the Schrodinger Opera- 
tors. Let us point out that standard Quantum Scattering Matrix S' 
on the line differs from our S by the multiplication on permutation 
matrix P from one side: S' = SP. For example, diagonal elements of S 
are equal to the Reflexion Coefficient, not to Transmission Coefficient as in 
standard matrix S'. The famous Reflect ionless Operators have Antidiagonal 
Scattering Matrix S in our sence for k = 2. For any k we immediately deduce 
from lagrangian property that 

<^t + Sjl^r^p + S pq^q >= ( S JP ~ S pj) < ^~ >= 

I 1 

< ?p + ,tp- >= V\ 2 -4 ^ 
Therefore our Scattering Matrix is symmetric. 

To prove unitarity we need to use reality of the Lagrangian plane. Con- 
sider now a real basis in the plane A: 

(J2tjiei)=T^ + )+TS(r) (29) 
i 

Its adjoint has a form 

TS(ip + ) +T(ip~) G A fc 
ip + = ip~ 

From that we conclude 

= 5 = 5* 

because we are coming to the same basis e G A. This line proves unitarity of 
S. Our Corollary is proved. 
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Remark 3 We can start with real basis in T\ of the form 



From the same arguments as above we deduce that A is a unitary matrix, and 



This is an imbedding of the space in Uk as a set of all symmetric unitary 
matrices. 

Let us describe now an Exceptional Discrete Spectrum. The following 
simple theorem is true. 

Theorem 4 Exceptional Discrete Eigenf unctions are completely defined by 
the eigenf unctions <p on the basis Y ' equal to zero in all nests (case of vertices). 



Here D 2 = D is the projection operator from functions on Y to the func- 
tions on Y' , putting all values outside subgraph equal to zero. 

Remark 4 In the case of edge operators we call 7 edge-nest 7 any edge R out- 
side ofY 7 touching some vertex-nest Pi. After that replacement all theorems 
containing the word 7 nest 7 remain true for edge operators. 

Proof of this theorem is very simple. Every such eigenfuction on the basis 
T' can be continuated to Y, taking zero value in all tails. This gives us an 
exact one-to-one correspondence between exceptional eigenvalues in Y and 
such special eigenvalues in Y', equal to zero in all nests. 

Corollary 3 The property of Operator L to have at least one exceptional 
eigenvalue has codimension not less than the number of nests or edge-nests 
(in the space of all real self adjoint Schrodinger Operators with finite domain 
where operator is different from the standard constant operator L ). In par- 
ticular, it is always not less than one, and Generic Operators have no excep- 
tional spectra. After generic small perturbation of L all discrete eigenvalues 
with |A| < 2 disappear. 



S = AA l e U k 




(30) 



0(P,) = O,Z = l,...,s 
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Proof of this corollary immediately follows from the fact that finite Graph 
T' has only a finite number of eigenvalues. Without the symmetry group, cor- 
responding eigenfunctions are generically nonzero in the nests (edge-nests). 

Remark 5 For the Graphs and operators with nontrivial symmetry group 
this simple counting of parameters does not work. Exceptional eigenvectors 
may necessarily appear in some symmetric cases-see examples 2 and 3 below. 

Remark 6 As Misha Gromov often explains in his lectures, Hyperbolic Ge- 
ometry is visible from the infinity as one- dimensional one. Therefore we may 
conclude, that for the discrete groups in 2D Lobachevski Plane with Noncom- 
pact Fundamental Domain of finite volume, Spectral Theory of the Laplace- 
Beltrami Operator should look in a sence 'similar' to the one on the graphs 
with k tails. We discussed this analogy with D.Kazdan, who pointed out to 
me that for arithmetic subgroups there are many discrete eigenvalues drown 
in the continuous spectrum. They disappear after nonarithmetic perturba- 
tion, as Peter Sarnak pointed out. In the case of graphs with k tails, we have 
simplified version of this picture for operators with symmetry: exceptional 
eigenvalues disappear after generic nonsymmetric perturbation. 

Example 2 Let k — 1 . Our Graph T is equal to finite subgraph T ', plus tail 
attached to it at one point (nest Pi). For generic graphs T, we have only 
Lagrangian plane T\, which is one dimensional, and no exceptional eigenval- 
ues. In the area where the operator is free (i.e. in the tail n > no), we have 
one real solution = aip + + bi])~ for every real A. In the Scattering Zone 
|A| < 2, we have 

= 0, -0 + = ip~ , b = a (31) 
|6/a| = l,Ae[-2,2] 

By definition of the Scattering Matrix, we have here 

s ( A) = b/a e Ui 

For the Zone of Normal Discrete Spectrum |A| > 2, we have an analytic 
continuation of scattering coefficient s(A) which has no sense of scattering 
anymore. Its poles a = give us points of normal discrete spectrum \ m . 
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// exceptional spectrum drown in the continuous one exists for this op- 
erator, we have an eigenfunction <fi' equal to zero in the tail. So we have 
a two-dimensional eigenspace and can consider Wronskian as a skew scalar 
product. 

Wronskian W((f>, 4>') for this exceptional value of A is equal to 
some finite cycle z G Hi(T'). For the graphs and operators with 
Z 2 -symmetry this possibility is generic, i.e., in general it cannot be 
destroyed by perturbation. 

To illustrate the last statement, consider Graph T with one tail and tri- 
angle T' = OAB attached to the vertex 0. Its vertices are numerated by 
nonnegative integers and letters A,B, P — n, . . . , 0, A, B with edges R n = 
[n, n — 1], Ra = [OA], Rb = [OB], Rab = AB. Take coefficients of the vertex 
Schrodinger Operator L in the form 

v n = 0, n > 0, v = u, v A = v,v B = w, 6 n:n _i = 1, n > 0, b 0:A = a, b 0]B = b, b A:B = c 

We have following set of equations: 

I u — A a b 
a v — A c 
\ c b w — A 

For the exceptional eigenvector we need such solution that ip± — ip — 0. It 
leads to the condition: 

X ex = w — bca^ 1 = v — acb~ l 

If \X ex \ < 2, we have an exceptional eigenvalue drown in the continuous 
spectrum. 

Z 2 -symmetry leads to the following coefficients: 

v = w,a = b 

Here we have \ ex = c. 





Example 3 Let us consider vertex operators for k = 2. We have two tails 
here and one or two nests in the Graph V . For almost all A we have here well- 
defined Monodromy Matrix as on the line. However, on the line monodromy 
was well-defined for all A. 
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Here we may have isolated "singular" values of A for which Lagrangian 
plane T\ is special: a monodromy map from the free basis of one tail into free 
basis of the another tail does not exist. It happens if our Lagrangian plane 
has a basis fa, fa such that fa = in the second tail and fa = in the first 
tail. 

Wronskian W(fa,fa) is equal to some finite cycle in V for such 
'singular' A = A*. 

It might happen even in the Scattering Zone as a generic possibility. Such 
singular values A = A* may be real, or they may appear by the complex adjoint 
pairs. 

To illustrate the last statement, we consider two cases: 
IT' is a triangle [OAB] where the vertex is a nest for both tails. 
2T' is a triangle [0i0 2 A] where the vertices 0i,0 2 are the nests. 
In the Case 1 we have vertices 

P = ...,n 1 ,...,0 1 =0,...,n 2 ,...,0 2 = 0,A,B 

and coefficients 

Kom-i = 1, n i > 0, i = 1, 2, b . A = a, b 0:B = b, b A:B = c, v = w, v A = u, v B = v 

As elementary calculation shows, we have solution fa equal to zero in the tail 
i = 1 or i = 2 iff 

(u-\*)(v-\*) = c 2 

Here we have real A* only, and W(ipi,ip2) = 
In the Case 2 we have vertices 

P = m, . . . , 0i, . . . , n 2 , ■ ■ ■ , 2 , A, m > 

and coefficients 

Kv.m-i = b n2 :n 2 ~i = l,rii> 0, &oi:o 2 = ^Mv-A = b,b 02 , A = c,u = v Ql ,v = v 02 ,w 

After elementary calculations we have 

\* = bca' 1 , W(fa, fa) = 6ca _1 [0i0 2 A] G H 2 (T, R) 

Here fa — in the tail with number i,i — 1,2, and fa{A) = 1. Exceptional 
eigenvalues with eigenf unctions equal to zero in both tails do not exist for this 
simple graph. 
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For the cases k > 3 the situation is much more complicated. 
Let us point out here that a very simple argument leads to the existence of 
discrete spectrum \ q > 4 for Laplace-Beltrami operators. Apply operator L 
to the set of ^-functions. Consider the square of norm 

(L^,L^) = \\L^\\ 2 

and the maximum of it for all ^-functions. We denote this maximum by Ml- 
For the Laplace-Beltrami operators, we shall use L = A + 2 as before. This 
choice of constant is optimal for this estimate. If this quantity is bigger 
than 4, we conclude that discrete spectrum exists for L. Moreover, 
this estimate is certainly nonexact. So, discrete eigenvalue \ q > 2 exists 
also if our estimate is exactly equal to 4 on the set of (^-functions. We are 
coming finally to the following estimates sufficient for the existence of discrete 
spectrum in most cases: 

1. Vert ices 

M L = max P (J2 b 2 P:P/ + V%) > 4 
p> 

For the vertex Laplace-Beltrami, we have Ml = maxp{mp + (mp — 2) 2 } 

2. Edges 

M L = max R (J2d 2 R .. R/ + Vl) >4 

R' 

For the edge Laplace-Beltrami, we have Ml = maxR{mp — 1 + mp/ — 1} 
where P U P' — R. It is easy to improve last inequalities, replacing the 
Operator L by L 2 . Sasha Veselov obtained better estimates for the case of 
vertices, for example, as he privately informed me. 

Factorization of Schrodinger Operators on Graphs 

Definition 7 We call Schrodinger Operator L acting on vertices or edges 
factorizable if it can be represented in the form 

L + C = QQ + ,C = const (32) 
L + C = Q + Q + U R 
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Factorization is special if function U R is equal to constant. For the case 
of vertices we consider only special factorization. Here operators Q, Q + are 
real adjoint to each other. An Operator Q+ maps functions of vertices into 
functions of edges by the following formula 

(Q + ^)r = T, c R--P^p ( 33 ) 
p 

(Q^)p = T, c R-P ( t ) R (34) 

R 

We call factorization formal if coefficients of operators Q, Q + are not real 
(and they are not adjoint actually to each other). 

Elementary substitution of this into our equations leads to the following 
result: 

Theorem 5 Representation of Operator L in the factorized form is equiva- 
lent to the set of equations: 

1. For vertices 

bp : p> = c R ..pc R ..p>, P ^ P' (35) 
PeR 

2. For edges 

d R -.R> = c R:P c R , :P , RDR' = P (36) 

V R + C= c\, p + c\, pl + Ur 

Corollary 4 For Operators L acting on edges: part of our factorization 
equations for the quantities c R: p through d R: p provides a complete set of al- 
gebraic equations for any given vertex P. This set of equations is overde- 
termined for such P that mp > 3. Compatibility conditions (in the 
form of algebraic constraint for the coefficients of operators) should 
be satisfied for factorization in this case. If all coefficients d R:P are pos- 
itive for all closed neighbors R ^ Rl ' , the solution c\. P is also positive. There 
is a formal factorization L = QQ f + Up such that the squares of coefficients 
c\. P are uniquely defined by the equation of factorization above. 
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Corollary 5 For operators acting on vertices: let all coefficients bp-p> are 
strictly positive for all closest neighbors P ^ P' . Take any finite contractible 
("tree-like") subgraph Y" in the Graph Y and its initial vertex (nest) P G Y" . 
Take any value of the constant C and any function Cr : p given along the 
boundary ofY" except the point Pq, where edges R look inside of our selected 
subgraph from the boundary vertices P. There is a special formal factorization 
of L with functions c 2 R . P uniquely defined by this data in the subgraph Y" . 
There is a value of the constant C depending on the subgraph Y" such that 
all quantities c R . P are positive. 

Proof of this easy follows from the form of the factorization equations: we 
solve the equations above for the quantities c R:P in both cases. Let me point 
out that factorization is purely local in the case of edges. So we proved our 
corollary for the case of edges. 

Consider the case of vertices more carefully. From the data on the bound- 
ary of subgraph Y" we can find unique real solution for the squares c R . P for 
any real constant C. However these quantities might take negative value. 
This leads to complex solutions in terms of CR.p. After that we find signs 
of Cr : p from the very simple part of our equation which does not contain 
squares, using initial data. 

To prove the last part of the Corollary we need to take constant C large 
enough and special initial data on the boundary of subgraph Y" . We take 
large enough values of c\. P on the boundary after choosing a large constant 
C . Here the edges R are attached to the boundary points P from inside of Y" . 
Solving the factorization equation in the direction to the endpoint Pq we shall 
go through the number of steps (layers) in Y", such that c R . P will be small, of 
the order for the edges looking from outside to the layer vertices. After 
that we shall find that the values of c Rl . P will be large and closed to C for 
the edges looking inside of layers. This ansatz is self consistent. Therefore we 
are coming to the desired positive solution. Corollary is proved. 

Conjecture: Let coefficients du-pr and their inverse d~£ R , for P ^ P', R ^ 
R' are positive and bounded. Let Graph Y is contractible and all numbers 
m P are also bounded. There is some positive value of C that Operator L 
admits a real factorization. 

Probably, this statement follows from the little inprovement of the same 
arguments as last Corollary. 
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Until now we did not investigate factorization for the graphs with non- 
trivial topology. 

Appendix: Two Remarks 
1. Nonlinear equations. 2. Fermionic Quadratic Forms 

We shall discuss here the continuous analog of our constructions for the 
case of vertices and nonlinear generalizations. After that we describe some 
properties of real fermionic quadratic forms. 

1. Consider any smooth manifold M with Riemannian metric and the 
linear selfadjoint Schrodinger Operator L acting in the space of the scalar 
functions. This operator can be obtained from the variational principle 

SW= I (J2ai k (x)^ + U(xM 2 )^)d n x 

J M 

in the Hilbert space of square integrable functions, where g(x) is determi- 
nant of the Riemann tensor in the local coordinate system x 1 , . . . , x n . What 
Gelfand told me in 1971 is that the expression 

/ {{14)<f> - {I4)$)yffa)d«x 

for the arbitrary pair of functions ip, is nonzero for the domain B. Using 
the Stokes formula, it can be reduced to the boundary dB = D. It leads to 
the integral 

j B ((L^<p-(L ( f ) )^^x)d n x = ^(0,^7 

Here d'-f means corresponding area element on the boundary. As we can see, 
what we get along the boundary is in fact an analog of Wronskian in our 
terminology. For the selfadjoint extension of operator L in the domain B 
you need to take boundary conditions in the form of some Lagrangian Plane 
in the last space of pairs (ip, 0) e H of functions on the boundary. 

In the main part of this work, we defined and used a discrete analog (for 
the Graph T) of the quantity W which is now the density of vector field 
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W j (x)y g(x)d n x on the manifold M in the continuous case. Our scheme 
corresponds in this case to the following: we take any 1-form Ui(x)dx l with 
compact support on M; W can be considered as a 1-chain on M ("current") 



For the pair of solutions (ip, <fi) of the equation Ltp = we are coming to 
the 1-cycle associated with current W. 

We did not clarify any continuous analogs of the cases k > 0. 

After the very useful discussion of our first work || with A.Schwarz in 
Maryland in December 1997, we came to the conclusion that this con- 
struction has a natural generalization to the nonlinear case: we 
may construct a closed H-valued 2-form on the space of solutions 
for nonlinear variational problems in both cases— continuous and 
discrete, i.e., on M and on the Graphs T. Here H = Hi pen (T, Z) for 
graphs, and H is a divergionless current for manifolds. Construction 
of this 2-form immediately follows from this work, applying it to the operator 
of second variation along the solution /, to the pair of variations for A = 0: 



This form is closed for A = 0. I will publish details in the next work. 

Remark 7 In one-dimensional case there is a natural class of continuous 
operators-the Schrodinger Operators on Graphs such that in every edge Ri a 
self-adjoint second order standard Schrodinger Operator Li is given, charac- 
terized by the set of real potentials Vi(x),x G Ri given in all edges R as the 
functions continuous in the closed edge (i.e. including boundary). For any 
vertex P we have k = nip edges Ri ending in it. Consider Symplectic Space 
Rp k given as a direct sum of 2-spaces R 2 associated with every edge Ri 



with canonical coordinates Pi,qi G R 2 . For any solution L^iix) = \ipi(x) 
in the edge R i} we have boundary values ip{P)^ip'{P) G Rf. Let for any 
vertex P a Lagrangian Plane Ap C Rp k is given. We call set of solutions 
Liipi(x) = Xtpi(x),x G Ri solution of the real selfadjoint operator L 




W f {8$,8<t>) = n 



Rf = R 2 © ... © R 2 
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on the graph iff in any vertex P full set of the boundary values 
belongs to the Lagrangian Planes Ap. 

We may replace operators L R by the Monodromy Matrices M R e SL 2 (R): 

M R : Rp — > Rp, , OR = P U P' 

All set of vertices can be considered as some kind of 'boundary' for the 
Schrodinger Operator L = UL R (union along all edges). The whole set of 
Lagrangian Planes A = UAp is some sort of 'self adjoined extension' for the 
Operator L. As P.Kuchment informed me in May 1998, these operators have 
been considered by mathematical and theoretical physicists for some concrete 
problems of Solid State Physics and Superconductivity in the late 80-s, who 
formulated self-adjoint junction condition in vertices in in the classical von 
Neumann-Krein form. The most popular conditions are following: 

l.^(P) = ^(P) = ^(P),$>W = a^(P) 

i 

24[{P) = ^(P) = ^'(P),J2MP) = <*j>'(P) 

Here indices i,j as before numerate all edges coming to the vertex P. (see 
B.Pavlov with N.Gerasimenko obtained also some results in the Scatter- 
ing Theory for such operators (see jTIj/j. We shall compare their results with 
our ideas in the next work. In the works of P.Kuchment and A.Figotin a 
pseudodifferential model on graphs was developed for photonic crystalls (see 
for example IfflJ). Certainly, no one of them discussed such problems in 
terms of Symplectic Geometry. 

Consider now any subgraph Y x C V such that dR C r x implies R C 
T. Therefore the subgraph r\ attached to the other part of the Graph Y by 
some edges R\,...,R P ending in the vertices Pi G Ti and in the vertices 
P/ G (TminusV) . For any self adjoint real Schrodinger Operator L on V we 
can define naturally the X-dependent Lagrangian Plane Apj in the Symplectic 
Space R 2p = R 2 Pl © ... © R 2 P generated by the values of the solutions = 
Xipi along the edges Ri and their first derivatives in the vertices Pi. All 
interaction of this subgraph with other part depends on this Lagrangian Plane. 
In a sence we may consider such subgraph equipped by the Lagrangian Plane 
as some more complicated kind of Selfinteracting Vertex with X-dependent 
boundary conditions. 



22 



2. Fermionic Quadratic Forms. Let us consider now a finite dimensional 
real linear space R n with basis and total space of external powers 

k=n 

A*(iT) = J2 Ah (R n ) 

Let us associate with every basic vector Cj a Fermionic Creation Operator 
ctj] we associate a vacuum vector 77 with unity 1 G i? = A C A*(i? n ): 

ej = 0,(77) 

The total space of exterior powers has a natural basis 

a ji ■ ■ ■ a jk 

(77), ji < ... <j k 

for all j,k. 

Definition 8 A Real Fermionic Quadratic Form is a selfadjoint operator L 
with real coefficients acting on the total space of exterior powers, written in 
the form: 

L = A pq a p a q + B pq a* p a q + C pq a*a* + const 

Here the Annihilation Operators a* are adjoint to the Creation Operators 
a p . They satisfy to the "canonical" relations 

CLpCLq (XqCLpj (XpClq \ CLqCLp PO. 

We have obviously 

A — — A R — R C - — A 

^pq ^qp-i J-'pq ^qpi ^pq PI 



Theorem 6 Take the new (noncanonical) basis af = a* ± aj An Operator 
L has the form 

L = D pq a p a q + const, D = A + B 
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Definition 9 A Bogolyubov Transformation or Canonical Transformation 
is an isomorphism of this Clifford Algebra given by the change of canonical 
basis: 

a P = Ppqbq + ^ Qpqbq 
a l = J2 QlP b 1 + E P ^ h *q 

Here new basis b*,bj satisfies to the same algabraic relations as the original 
one (both of them should be canonical). 



Theorem 7 For any Canonical Transformation both matrices 0± = P ±Q 
are orthogonal 0± G O n (R). The Transformation Rule for the matrix D is 
following: 

D = O.D'O- 



where 



L = D pq a+a q = D' rs b+b~ s 



This theorem can be easily proved by the direct elementary verification. As 
a corollary from this we have 

Corollary 6 An operator L can be reduced to the diagonal form using Bo- 
golyubov Transformations above and Transformation Rule for matrix D = 
A + B. Eigenvalues fij 1 j 2 ...j k of the Operator L in the space A*(R n ) can be 
computed in the following way: 

Mil + ---H k I'.i—j- ■ < ••• <Jk 
where Hj are the eigenvalues of the 'Absolute Value Operator' 



\L\ = VI*L 

in the space R n 

This algebraic statement has been found as a lemma in the present au- 



thor's Appendix to our joint work with Misha Shubin (see [[T0[)- 111 this 
Appendix I developed a nice analytic Witten-like approach trying to find 
right analog of the Morse Inequalities for the Vector Fields on real mani- 
folds: necessity to diagonalize an arbitrary Real Fermionic Quadratic Form 
appears naturally here. 
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(Let us remind that in the Witten's approach to the Morse inequalities 
for ordinary functions we have A pq = 0; therefore only one orthogonal matrix 
has been used for the diagonalization). Several experts in Quantum Field 
Theory pointed out to the present author that this elementary observation 
has never appeared in the literature, so I decided to repeat it here once more. 
(The last one was A.Schwarz who told me that in the December 1997.) 
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